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A magnetic field applied to type-II superconductors introduces quantized vortices that locally 
quench superconductivity, providing a unique opportunity to investigate electronic orders that may 
compete with superconductivity. This is especially true in cuprate superconductors in which mu¬ 
tual relationships among superconductivity, pseudogap, and broken-spatial-symmetry states have 
attracted much attention. Here we observe energy and momentum dependent bipartite electronic 
superstructures in the vortex core of BLS^CaCi^Os+s using spectroscopic-imaging scanning tun¬ 
neling microscopy (SI-STM). In the low-energy range where the nodal Bogoliubov quasiparticles are 
well-defined, we show that the quasiparticle scattering off vortices generates the electronic super¬ 
structure known as “vortex checkerboard”. In the high-energy region where the pseudogap develops, 
vortices amplify the broken-spatial-symmetry patterns that preexist in zero field. These data reveal 
canonical d- wave superconductivity near the node, yet competition between superconductivity and 
broken-spatial-symmetry states near the antinode. 


The electronic states in cuprates exhibit distinct fea¬ 
tures depending on energy and momentum as schemati¬ 
cally shown in Fig. 1A [1, 2]. The low-energy near-nodal 
states host the homogeneous d -wave superconductivity 
that manifests itself in the Bogoliubov quasiparticle inter¬ 
ference (BQPI) patterns imaged by SI-STM [2-5]. BQPI 
is no longer observed above the doping-dependent extinc¬ 
tion energy Ao and outside the diagonal line in momen¬ 
tum space connecting (7r/ao,0) and (0, 7r/ao), where ao 
denotes Cu-O-Cu distance [2]. The higher energy states 
near the antinode are governed by the pseudogap whose 
apparent magnitude Ai is spatially inhomogeneous and 
the quasiparticle excitations near Ai break rotational 
and translational symmetries of the Cu0 2 plane [2, 6, 7]. 

In order to establish the relationship among these elec¬ 
tronic states, it is indispensable to investigate how the 
pseudogap and the broken-spatial-symmetry state are af¬ 
fected when the superconductivity is suppressed. Intro¬ 
duction of vortices is one of the ways to suppress su¬ 
perconductivity. It has been shown in La- and Y-based 
cuprates that the electronic orders that break the spa¬ 
tial symmetry are enhanced or even generated by vor¬ 
tices [8-13]. However, the detailed energy, spatial and 
momentum-space structures of the vortex-enhanced or¬ 
ders are unknown. To address this issue, we utilize SI- 
STM in the vortex state owing to the following three ad¬ 
vantages. First, vortices can suppress superconductivity 
at the lowest temperatures where the thermal broaden¬ 
ing effect is negligible, making it possible to study the 
precise energy scale of the phenomenon. Second, atomic- 
scale spatial resolution of SI-STM is highly beneficial not 
only to identify the locations of vortices but also to de¬ 
termine the real-space structure of vortex-induced states. 
Finally, by using the Fourier transformation, SI-STM ac¬ 
quires the momentum-space resolution even under mag- 
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netic fields that enables us to discuss the near-nodal and 
antinodal states separately. 

We choose optimally doped Bi 2 Sr 2 CaCu208+<5 (super¬ 
conducting transition temperature T c ~ 90 K) as a sam¬ 
ple because of its high quality surface necessary for SI- 
STM [14]. Pioneering SI-STM studies of the vortices in 
Bi-based cuprates have discovered that an electronic su¬ 
perstructure, so-called “vortex checkerboard”, is nucle¬ 
ated in the vortex core [15-17]. Although a possible 
connection between the “vortex checkerboard” and the 
electronic order that breaks spatial symmetry has been 
discussed [15-17], the electronic states in the vortex core 
are still elusive, probably because the energy ranges so 
far studied are mostly below Aq and little is known about 
the momentum-space electronic states in the vortex state. 
Here we reveal energy-dependent bipartite nature of the 
electronic superstructures in the vortex core using the 
Fourier-transform SI-STM over a wide energy range. 

Figure 1, C to F show differential tunneling conduc¬ 
tance g( r, E, B ) maps at energy E = ±10 meV taken in 
magnetic fields of B = 0 T and 11 T in exactly the same 
field of view. Here, r denotes the position on the sur¬ 
face. Vortex cores are identified as enhanced g{r,E,B) 
regions with the “vortex checkerboard” structure [15-17]. 
As shown in Fig. IB, the vortex alters the spectrum in 
two different energy regions [16-20]: the emergence of 
conductance humps around \E\ ~ 10 meV < A 0 and 
the suppression of the peaks at Ai. Since these two en¬ 
ergy regions occupy different sectors in momentum space, 
near-nodal and antinodal states (Fig. 1A), it is intriguing 
to explore the momentum-space characters of the vortex- 
induced electronic states. 

By taking the Fourier transformation from the spec¬ 
troscopic images, we can estimate the characteris¬ 
tic wavevectors q(E,T>)’s of electronic-state modula¬ 
tions. However, in heterogeneous systems such as 
Bi 2 Sr 2 CaCu 2 0 8 +<5, g(r,E,B) not only reflects the r de¬ 
pendence of the local density-of-states (LDOS) at E 
but also includes LDOS modulations at different ener¬ 
gies because of the r-dependent tip elevation associated 
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FIG. 1: Nodal-antinodal dichotomy and spectroscopic features of vortices. (A) Schematic illustration of the excitation 
gap in momentum space showing the dichotomy between the d -wave superconductivity near the node (light blue area) and the 
antinodal states governed by the pseudogap (yellow area). These two regimes are separated by the line connecting (7r/ao,0) 
and (0,7r/ao). (B) Comparison between tunneling spectra taken at B = 0 T (open symbols) and 11 T (solid symbols). The 
measurement temperature is 4.6 K. Red and blue data depict the spectra spatially averaged over the regions near vortices 
and far from vortices, respectively. Detailed point spectra before averaging are shown in fig. S3. (C and D) Differential 
conductance maps g(r,E,B) at E — +10 meV in magnetic fields B = 0 T and 11 T, respectively. White arrows in (C) 
denote the Cu-0 boding directions. Vortices and their internal structures (“vortex checkerboard”) are clearly imaged in (D). 
(E and F) Differential conductance maps g(r,E,B) at E — —10 meV in magnetic fields B = 0 T and 11 T, respectively. The 
tunneling conductance at each location was obtained by numerical differentiation of the current-voltage characteristics and by 
post-smoothing with the energy window of ±2 meV. 


with the feed-back loop. This so-called set-point ef¬ 
fect can be suppressed by taking a ratio Z(r,E,B) = 
g(r, +\E\, B)/g(r f — \E\, B), which faithfully represents 
the ratio of the LDOS at E\E\ [21, 22]. 

First, we focus on the low-energy near-nodal region 
and argue the origin of the “vortex checkerboard”. At 
B = 0T, the only relevant phenomenon near the node is 
the BQPI that is described by the “octet model” [4, 5] 
in which the eight tips of the banana-shaped constant- 
energy contours in momentum space dominate the quasi¬ 
particle scatterings, resulting in a set of energy-dispersive 
characteristic wavevectors (i = 1 , 2,- • •, 7)(See sec¬ 
tion 2.1 of supplementary online text). Figure 2A de¬ 
picts the typical BQPI pattern seen in Z q (c\,E,B = 
0 T), the Fourier-transformed image of Z(r,E,B = 
0 T), showing the octet q^’s. (In the optimally doped 
Bi 2 Sr 2 CaCu 208 ++ signals at q 4 and qs are weak in 
Z [2, 23].) We have performed standard BQPI analy¬ 
sis [2, 5] (See section 2.1 of supplementary online text), 
and confirm that the BQPI is restricted in the near-nodal 
region below the extinction energy Aq ^ 30 meV [2]. 


An important question here is whether the “vortex 
checkerboard”, which is most prominent around |^7| ^ 
10 meV < Aq, represents an electronic order or not. 
To answer this question, we have repeated the same 
Z(y,E,B) analysis at B = 11 T in the same field 
of view (See section 2.1 of supplementary online text). 
As shown in Fig. 2B, no additional peak is detected 
whereas the intensity of each q^ peak depends on the 
magnetic field. Figure 2C highlights the field-induced 
change obtained by subtracting Z q (q,E,B = 0 T) from 
Z q (q,E,B = 11 T). The enhanced intensity appears 
at qi, which represents the wavevector of the “vortex 
checkerboard”. We note that the intensities at q 2 , q 3 , q6, 
and q 7 are suppressed by the magnetic field. All of these 
scattering q^’s reverse the sign of the d-wave supercon¬ 
ducting gap between the initial state and the final state, 
while the sign is preserved in the the case of qi scattering. 
Such suppression and enhancement of sign-reversing and 
sign-preserving scatterings, respectively, are exactly what 
are expected from the coherence factors of the quasi¬ 
particle scatterings off vortices [24-26], suggesting that 
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FIG. 2: Magnetic-field effect on the low-energy state 
showing the BQPI origin of the “vortex checker¬ 
board”. (A and B) Two-fold symmetrized Fourier- 
transformed images Z q ( q, E , B ) taken in the field of view of 
470 x 470 A 2 . The tunneling conductance at each location 
was taken by the standard lock-in technique with the modu¬ 
lation amplitude of 2.5 mV rms . The measurement tempera¬ 
ture is 4.6 K. Red and blue arrows indicate the sign-reversing 
(q 2 , q3, q 6 , qr) and sign-preserving scattering (qi, q 4 , qs) 
wavevectors, respectively. Black circles show the Bragg spots. 
(C) Difference between (A) and (B); Z q { q, E — 15 meV, B = 
11 T) — Z g (q, 5 = 15 meV, B = 0 T). (D and E) Energy- 
dependent line profiles taken along the black arrows in (A) 
and (B), respectively. Red open circles denote the position of 
qi peak that has been determined by fitting the line profile 
at each energy by Lorentzian function. Error bars indicate 
the full-width-half-maximum of the qi peak. Dashed red line 
in (D) and solid red line in (E) are the dispersions calcu¬ 
lated from q2, q3, q 6 , q7 based on the “octet model”. Note 
that the vortex-induced signal exhibits the dispersion that is 
consistent with the BQPI. (F) Energy-dependent line profile 
taken along the black arrow in (C), showing the magnetic- 
field induced change. See section 2.1 of supplementary online 
text for details. 


the vortex-enhanced BQPI is the cause of the “vortex 
checkerboard”. 

The BQPI scenario is further supported by the energy 
dispersion in the vortex-enhanced signal at q 4 (Fig. 2, 
D and E). The observed qi dispersion agrees well with 
the behavior calculated from q 2 , q 3 , q6, and q 7 based 
on the octet model. Together with the fact that the 
vortex-enhanced signal diminishes near Aq (Fig. 2F), we 
ascribe the “vortex checkerboard” to the BQPI and con¬ 


clude that no electronic order is nucleated in the vortex 
cores at E < Aq. We note that other Friedel-type oscil¬ 
lations such as bound-state oscillations in the quantum- 
limit vortex core [17] may also be relevant (See section 
2.2 of supplementary online text). 

Next, we study the effect of vortices near A 4 . The elec¬ 
tronic feature here is characterized by the bond-centered 
unidirectional electronic entity that breaks both rota¬ 
tional and translational symmetries [2, 6, 7, 21], being 
reminiscent of the short-range charge order detected by 
X-ray scattering [27, 28]. For the purpose of brevity, we 
call this electronic entity as “nanostripe” hereafter. To 
visualize the “nanostripe” which is most prominent at 
Ai, we follow the procedure used in [2, 6]. Since Ai is 
spatially inhomogeneous, we first normalize E by local 
Ai(r) and map Z( r,e = F/Ai(r) = 1,5). Figure 3A 
shows Z(r, e = 1, B = 0 T) in the same field of view as 
Fig. 1. The “nanostripe” in the optimally-doped sam¬ 
ple is weak in intensity [23] and is observed only in the 
limited regions where Z( r, e = 1,5) is larger. We find 
that these regions have larger Ai(r) and vortices tend 
to reside there, suggesting a vortex pinning mechanism 
associated with the pseudogap (See section 2.3 of supple¬ 
mentary online text). We have repeated the same mea¬ 
surement at B = 11 T (Fig. 3B) and have revealed that 
Z( r, e = 1,5) is enhanced in the vortex cores (Fig. 3C). 
As shown in the insets of Fig. 3, A to C, the structure 
of the “nanostripe” that preexists in the absence of the 
vortex core is unchanged but its contrast is enhanced. 

It has been known that the “nanostripe” consists of 
two sets of wavevectors [6]: Q X:V = (27r/ao,0), (0,27r/ao) 
whose inequivalent intensities represents the degree 
of broken rotational symmetry and y ~ (3/4 x 
(27r/a 0 ),0),(0,3/4 x (27r/ao)) that features the broken 
translational symmetry. To test which broken symme¬ 
try is affected in the vortex core, we perform the Fourier 
analysis. By applying a mask generated from the image 
of vortices, we restrict our field of view in the vicinity of 
the vortex core to effectively extract the vortex-enhanced 
features [14]. As shown in Fig. 3, D and E, Fourier peaks 
corresponding to Q and S are identified at both 5 = 0 T 
and 11 T. In the difference image shown in Fig. 3F, 
the intensity at S is enhanced, whereas the change at 
Q is small. This observation indicates that the vortex 
core predominantly amplifies the translational-symmetry 
breaking. 

The bipartite vortex-induced changes observed at E < 
A 0 and at E ~ Ai provide important implications to 
understand the electronic-state competition in cuprates. 
First, the BQPI origin of the “vortex checkerboard” at 
E < Ao suggests that canonical phenomenology of d- 
wave superconductivity applies in the near-nodal region 
even when the vortices are introduced. 

At 5 ^ Ai near the antinode, the “nanostripe” 
amplified in the vortex core resembles the magnetic- 
field-enhanced charge orders in Y-based cuprates [9-13]. 
This indicates that the magnetic-field-enhanced broken- 
spatial-symmetry state is universal in cuprates. 
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FIG. 3: Magnetic-field effect on the “nanostripe” at the pseudogap energy scale. (A and B) Conductance ratio 
maps Z( r, e = 1 , B) in magnetic fields B — 0 T and 11 T, respectively. These images are taken in the same field of view 
as Fig. 1, (C) to (F). The tunneling conductance at each location was obtained by numerical differentiation of the current- 
voltage characteristics and by post-smoothing with the energy window of ±10 meV. The measurement temperature is 4.6 K. 
(C) Difference between (A) and (B); Z(r, e — 1,B = 11 T) — Z(r, e = t,B = 0 T). Insets of (A) to (C) are the magnified 
images of the regions marked by boxes in the main figures. Red and blue lines in the insets indicate the directions of the 
“nanostripe”. Dots denote the locations of the Cu atoms. (D to E) Fourier-transformed images Z q ( q, e = 1 , B) in magnetic 
fields B = 0 T and 11 T, respectively. The field of views are restricted in the vicinity of the vortices by applying a mask [14]. 
(F) Difference between (D) and (E). In these Fourier-transformed images, dotted and solid circles indicate the locations of 
Q x ,y = (27r/ao, 0), (0, 27r/ao) and S XjV ^ 3/4 x (27r/ao, 0), 3/4 x (0, 27r/ao), respectively. 


It is noteworthy that even though the contrast of the 
“nanostripe” seen in Z(r, e = 1, B) is enhanced near the 
vortex core, the spectral weight at |£7| ~ Ai is suppressed 
as shown in Fig. IB. If the missing weight is caused by the 
reduction of electrons condensed in the pseudogap state, 
the enhanced “nanostripe” in the vortex core should be 
distinct from the pseudogap formation. This picture is 
unlikely because the previous SI-STM results in the ab¬ 
sence of vortices show close correlation between the pseu¬ 
dogap and the “nanostripe” [2, 23, 29]. Rather, we infer 
that Bogoliubov quasiparticles exist even near the antin¬ 
ode at |^7| ~ Ai and that the the missing weight is origi¬ 
nated from the suppression of superconductivity [17, 31]. 
Since the pseudogap or the “nanostripe” also develops 
at the same energy scale and two different orders do 
not share the same region in energy-momentum space, 
it is plausible that superconductivity and the pseudo¬ 


gap are nearly degenerated but competing states near 
the antinode [9, 13, 30, 31]. This scenario is consis¬ 
tent with the heterogeneous electronic states in real space 
because randomly-dispersed dopants and defects locally 
perturb the balance between superconductivity and the 
pseudogap, and consequently bring about the nano-scale 
mixture of them. The observed vortex-enhanced “nanos¬ 
tripe” suggests that the vortex core is one of these per¬ 
turbations that affects the competition near the antin¬ 
ode. Since the energy scales of superconductivity and 
the pseudogap near the antinode change with doping, we 
anticipate that the vortex-core electronic states at dif¬ 
ferent doping levels provide more insights into the rela¬ 
tionship among superconductivity, pseudogap, and the 
“nanostripe”, leading to establishing the phenomenology 
of cuprates. 
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Supplementary Materials 


unaffected. 


1. MATERIALS AND METHODS 
1. 1. Procedure of SI-STM experiments 

SI-STM experiments were performed at a tempera¬ 
ture of 4.6 K with a commercial low-temperature ultra- 
high-vacuum STM (Unisoku USM-1300) modified by our¬ 
selves [32]. Single crystals of Bi2Sr2CaCu208+<5 were 
grown by the traveling-solvent floating zone method and 
were annealed to have optimal hole concentration. The 
superconducting transition temperature T c ~ 90 K was 
determined by the magnetization measurement. The 
sample was cleaved in-situ at ^77 K to obtain a clean 
and flat (001) surface and was transferred quickly to the 
STM unit kept at 4.6 K. We used an electro-chemically 
etched tungsten wire as an STM tip, which was cleaned 
and characterized in-situ with a field-ion microscope. All 
the SI-STM data were taken with the feed-back set point 
at a sample bias voltage of -150 mV and a tunneling cur¬ 
rent of 150 pA. Magnetic field was applied perpendicular 
to the cleaved (001) surface. Whenever we changed the 
magnetic field, the sample was heated up to ^30 K to 
make the vortex distribution inside the sample uniform. 


1. 3. Procedure to generate masks 

Vortices are clearly seen in the difference map 
5g(r,E = +10 meV, D = 11 T) = g(r, E = 

+10 meV, B = 11 T) — g(r,E = +10 meV, B = 0 T) 
(Fig. S2A). We first apply low-pass filter to the 5g( r, E = 
+10 meV, B = 11 T) map with a cut-off wavelength of 
q ~ 0.05 x 27r/ao- Constant contours of this filtered image 
are used for the boundaries of the masks (Fig. S2B). 

The ‘vortex’ and ‘matrix’ regions used to examine the 
vortex-induced change are shown in Fig. S2C. The spec¬ 
tra shown in Fig. IB of the main text are the spectra 
averaged in these regions. In Fig. S3, we depict the de¬ 
tailed point spectra near the vortex before averaging. 

The mask used for the Fourier analyses (Fig. 3, D to 
F in the main text) and the restricted Z(r, F,F) maps 
at B = 0 T and 11 T are indicated in Fig. S2D, E, and 
F, respectively. 

2. SUPPLEMENTARY TEXT 
2. 1. BQPI analyses based on the “octet model” 


1. 2. Procedure to obtain Z(r,E,B) and Z(r, e, B) 

It is well-known that electronic states of 
Bi 2 Sr 2 CaCu 208+(5 are heterogeneous especially near 
Ai. Therefore, in order to argue the vortex-induced 
local changes in the electronic state, it is indispensable 
to compare two SI-STM data sets, with and without 
a magnetic field, in exactly the same field of view. 
This is challenging because the actual SI-STM images 
are inevitably distorted in an uncontrollable manner 
due to the creeping of the piezoelectric scanner, etc. 
In order to obtain distortion-free images from the 
observed ones, we utilize the so-called Lawler-Fujita 
algorithm [6], which estimates the local distortion as a 
phase shift in the crystal-lattice modulations. We first 
estimate and correct the distortions in the topographic 
images simultaneously taken with the g(r, F,F) images 
and correct the spectroscopic images, g( r, F, B ) and 
F(r,F, F), using the same local distortions. Various 
images obtained by SI-STM after the correction are 
shown in Fig. SI. 

In order to obtain Z( r, e, F), we have to estimate Ai(r) 
at each pixel. We define Ai(r) at the energy of the 
peak in the unoccupied state in the g(r, F, F = 0 T) 
spectrum and normalize the energy as e = E/ Ai(r). If 
there are multiple peaks in the spectrum, we have cho¬ 
sen the highest peak. We have used Ai(r) estimated 
at B = 0 T to normalize the energy in a field as well. 
Even if we use Ai(r) estimated at B = 11 T to calcu¬ 
late Z(r, e, B = 11 T), the conclusion in this paper is 


BQPI brings about LDOS modulations which would be 
reflected on g( r, F, B). However, in actual fact g( r, F, B ) 
may also contain extrinsic modulations caused by the set- 
point effect [21, 22]. The set-point effect can be largely 
suppressed in F(r,F,F), which we utilize throughout 
this work. The analysis of Z(r, F,F) has another ad¬ 
vantage in that Z(r,F,F) picks up the particular mod¬ 
ulations in which q(+|F|) = q( —|F|) and the phase dif¬ 
ference between the q(+|F|) and q(—|F|) modulations 
is large (~ i r). These features are exactly what are ex¬ 
pected in the BQPI [22, 33]. 

We use the “octet model” to analyze the BQPI seen 
in Z g (q, F,F); the Fourier transformed F(r,F,F) [3- 
5]. Here, the seven scattering vectors q^ (i = 1, 2,- • •, 
7) connecting the eight tips of the banana-shaped con¬ 
stant contours in momentum space govern the BQPI be¬ 
cause the joint density of states takes maximum for these 
wavevectors (Fig. S4A). As shown in Fig. 2, A and B of 
the main text and in Fig. S5, a set of energy-dispersive 
wavevectors are detected in Z q (q, F, B) and each of the 
wavevectors can be assigned to one of the octet wavevec¬ 
tors q^. The intensities at q 4 and qs are weak, being 
consistent with the previous reports [2, 23]. 

To determine the precise locations of q^’s at B = 0 T 
and 11 T, we fit the peaks in Z g (q, F, B) with the two- 
dimensional Lorentzian plus linear background; 


f(Qx,q y ) = fo+ 


^ '!■>: QxO ^ | QyO ^ _j_ 


■ Ec x q x -\-c y q y , 


(SI) 
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where /o, c x , and c y are the fitting parameters associated 
with the linear background, whereas A, q x o, and 

G y are the fitting parameters correspond to the ampli¬ 
tude, the q x and q y components of the peak location, 
and the q x and q y components of the half-width-at-half- 
maximum of the peak, respectively. Figure S4B shows 
the energy dependence of the absolute value of the ob¬ 
served q^’s. The signals at q 2 , q 3 , q6, and diminishes 
at about 30 meV, which set the extinction energy A 0 
(Fig. S4, B and D). At E > Ao, signals at qi and qs are 
still there but lose their energy dependence. The signal 
near qs turns into S, which is one of the ingredients of 
the “nanostripe”. 

We obtain the normal-state Fermi surface and the su¬ 
perconducting gap dispersion using pairs of (q 2 , qe) and 
(q 3 , < 17 ) • Figure. S4C depicts the Fermi-surface loci that 
sustain coherent Bogoliubov quasiparticles. They are 
limited inside the diagonal line connecting (7r/ao,0) and 
(0, 7r/ao) [2] and are hardly affected by a magnetic field. 
By contrast a magnetic field suppresses the near-nodal 
superconducting gap as shown in Fig. S4D. These fea¬ 
tures are consistent with the behaviors observed in the 
different cuprate superconductor Ca2- ;r Na ;r Cu02Cl2 [24] 
and can be associated with the Volovik effect [34]. 

In order to check the validity of the “octet model” 
especially in a magnetic field, we have performed the fol¬ 
lowing analyses. First, we fit the Fermi-surface loci to 
the quarter circle. Next, the superconducting gap dis¬ 
persions A (6) are fitted by the following d-wave form 
with an higher-order term. 

A (9) = A B qpi[7Lcos(2$) + (1 — A) cos(60)]. (S2) 

Here, 9 represents the Fermi-surface angle defined in 
Fig. S4C and A B qpi and A are fitting parameters. Using 
the obtained fitting parameters, the energy dispersions 
of all the qds can be calculated and plotted in Fig. S4B 
and in Fig. 2, D and E of the main text. The “vor¬ 
tex checkerboard” is characterized by the field-enhanced 
signal at qi. We note that the energy dispersion of qi 
at B = 11 T well coincides with the calculated one, as 
in the case at B = 0 T. Together with the fact that the 
q 4 and qs are also consistent with the calculated disper¬ 
sions, we conclude that the field-induced change in the 
electronic state can be explained in the framework of the 
“octet model” as long as E < Aq. 


2. 2. Possible origins of the enhanced qi 
modulations 

The applicability of the “octet model” implies that the 
observed electronic-state modulations are associated with 
the Bogoliubov quasiparticles that reside on the near- 
nodal Fermi surface. There is more than one origin to 
cause such Friedel-type oscillations. Most naively, the en¬ 
hanced quasiparticle scattering off vortices may result in 


the enhanced qi [24]. Another possible origin is the spa¬ 
tial oscillations of the vortex bound state in the quantum- 
limit vortex core [35], as proposed by Yoshizawa and 
coworkers [17]. This model naturally explains the peaks 
in the spectrum at ~ ±10 meV as the discrete bound 
states. Nevertheless, more studies are necessary to verify 
the validity of the bound-state scenario, because it is not 
clear whether the vortex of a cuprate is in the quantum 
limit and the formation of the vortex bound state in a d- 
wave superconductor is still a controversial issue [36, 37]. 

2. 3. Vortex pinning and the competition between 
the “nanostripe” and superconductivity 

It is interesting to examine the relationship between 
the locations of vortices and the various electronic het¬ 
erogeneities, since this comparison may give us a hint 
to identify the elementary process of the vortex pinning. 
Although previous SI-STM studies have suggested that 
vortices in BUS^CaC^Og+tS tend to be pinned at the re¬ 
gions where the Ai(r) is large [17, 38], these experiments 
were conducted only in a field and thereby could not ex¬ 
clude the possibility that the the pseudogap itself would 
be influenced by vortices. Here we compare the locations 
of vortices with the zero-field electronic heterogeneity at 
Ai(r). 

Using the same procedure described in supplementary 
section 1.3, we first identify the locations of vortices 
(Fig. S6, A and B). We examine the azimuthally aver¬ 
aged cross-correlation function (Fig. S6C) between the 
smoothed vortex map Figs S6B and three different spec¬ 
troscopic images at Ai(r): spatial variation of Ai(r) it¬ 
self (Fig. S6D), Z( r, e = 1, F>) (Fig. S6E), and the local 
amplitude of the broken-translational-symmetry state of 
the “nanostripe”, which is nothing but the local ampli¬ 
tude of the modulations at S x , y (Fig. S6F). Here, we 
define the local amplitude of these modulations As (r) as 
follows, 

-ls(r) = .l(S.,..r) I .l(S f/ . r). 

A(S v ,r) = y]Z(r')e* s "' r /A(r',-r) 

r' 

« ±=^^(S,-k)e ik ' r e- k2 / 2A2 , {v = xo t y), 

V N k 

(S3) 

where /a(U = (A 2 /27r)e — a 2 l r l 2 / 2 , 1/A is the cut-off 
length scale, and Z is the complex Fourier transform of 
Z(r, e = 1, B). As shown in Fig. S6C, all of these quanti¬ 
ties exhibit strong correlations with the locations of vor¬ 
tices. Since vortices are generally pinned at weakly su¬ 
perconducting regions, these results suggest that at least 
one of these quantities would represent the fundamental 
measure of the weakness of superconductivity. Although 
the microscopic mechanism of vortex pinning is unclear 



at present, the observed correlation between the vortex is prominent, indicating the competition between super¬ 
location and As(r) clearly indicates that the supercon- conductivity and the “nanostripe” in Bi2Sr2CaCu208+<5. 
ductivity is weak in the region where the “nanostripe” 



9 


3. SUPPLEMENTARY FIGURES 


Topography T(r, 8) g (r, E = +10 meV, 8) Z (r, E = 10 meV, 8) Z (r, e = 1,8) ^i(r, B) 



FIG. SI: Top and middle rows display two series of SI-STM results on a 380 x 380 A 2 field of view at B = 0 T and 11 T, 
respectively. Bottom row indicates difference images, which are obtained by subtracting the results at 0 T from those at 11 T. 
First column; topographic images, second; g(r,E = +10 meV, B), third; Z(r,E = 10 meV,B), forth; Z( r, e = 1, R), and fifth; 
Ai(r) . Although an apparent enhancement of Ai(r) is observed in some of the vortex cores, changes in the original g{ r, E, B ) 
spectra near Ai(r) are very small; due to small dg( r, E, B)/dE near Ai(r) , small change in g( r, E, B) gives rise to large effect 
on the estimation of Ai(r) . See, for example, middle spectra along the line C8 in Fig. S3. We have used Ai(r) at B — 0 T 
to create both of Z( r, e = 1,B = 0 T) and Z( r, e = 1,B = 11 T) maps. Black arrows in the topographic image at B = 0 T 
denote the Cu-0 bonding directions. 


































































10 


5g(r,E=+10 meV) LPF-5g(r,E = +10 meV) 


A 



e 



QlM 

• 


* 

0 M 

i I 

+0.4 


High 

f ^ 

■ 

■ ^ 

j£. ■ 

ia* 


-0.4 

50A 

Low 



Red : Vortex region 
Blue : Matrix region 




FIG. S2: (A) Difference between the conductance maps at B = 0 T and 11 T; Sg( r, E = +10 meV) = g( r, E = +10 meV, B = 
11 T) — g( r, E = +10 meV, B — 0 T). (B), Low-pass filtered image of (A). A cut-off wavevector is chosen to be 0.05 X 2ir/ao 
in the filtering process. Corresponding spatial resolution is indicated as a black circle and an arrow. Red dashed, blue dashed, 
and yellow lines are the contours with the value of 65, 8, 35% of the difference between the maximum and minimum values of 
this image, respectively. (C) Separation of “vortex” and “matrix” regions for the spatially averaging process used in Fig. IB 
of the main text. Red and blue regions denote “vortex” and “matrix” regions and are defined as the regions surrounded by 
the red and blue dashed lines in (B). (D) A mask for the effective extraction of the field-effect on the “nanostripe”. White 
and black regions denote the “vortex” and “matrix” regions, which are separated by the yellow contours in (B). (E) and (F) 
Masked conductance ratio maps at PG energy at B = 0 T and 11 T, respectively. Fourier transformed images of these masked 
maps are shown in Fig. 3, D and E in the main text, respectively. 
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g (r, E = +10 meV, 8 = 11 T) 
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FIG. S3: Continued on next page. 
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Bias voltage (mV) Bias voltage (mV) Bias voltage (mV) Bias voltage (mV) Bias voltage (mV) 

FIG. S3: Line profiles of tunneling spectra along 18 paths composed of 10 paths on “crests” and 8 paths on “troughs” of the 
“vortex checkerboard” seen in g(r,E = +10 meV, B = 11 T). Bottoms of these line profiles correspond to the starting points of 
the arrows drawn in the vortex image. Line profiles on the crests are shown in this page. Line profiles on the troughs are shown 
in the previous page. Blue and red curves in these profiles denote the tunneling spectra at B = 0 T and 11 T, respectively. 
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FIG. S4: (A) Schematic illustration of the momentum-space electronic structure. White and gray regions indicate the region 
with the opposite sign of d -wave superconducting gap. Black and gray solid lines represent normal-state Fermi surface and 
constant energy contours of the Bogoliubov quasiparticle dispersion. Green bold lines denote constant energy contours at a 
representative energy. Black dashed circles correspond to the eight high-density-of-states regions, which primarily contribute 
the BQPI at a given energy. Red and blue arrows indicate sign-reversing and sign-preserving scattering vectors [24]. (B) 
Energy dependence of the absolute value of the observed qPs at B = 0 T (open symbols) and 11 T (solid symbols). Dashed 
and solid lines denote the expected dispersions at B = 0 T and 11 T, respectively. (See Supplementary Section 2.1 for details.) 
(C) Fermi surface loci obtained from BQPI peak locations q i(E), at B = 0 T (open symbols) and 11 T (solid symbols). Green 
and orange circles in right hand side of this panel are obtained from pairs of (q 2 , qe) and (q 3 , qr), respectively. Symbols in 
left hand side represent the average of them. Blue dashed and red solid lines denote the results of quarter circle fitting. (D) 
Superconducting gap dispersion at B = 0 T (blue symbols) and 11 T (red symbols). Green and orange circles in right hand 
side of this panel are the dispersions estimated from pairs of (q 2 , qe) and (q 3 , qr), respectively. Symbols in left hand side 
represent the average of them. These are fitted to A(0) = Abqpi[Acos(20) + (1 — A) cos(60)] as shown by the solid curves. 
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FIG. S5: Top and middle rows display two series of Z q (q, E, B) at B = 0 T and 11 T, respectively. Bottom row indicates 
difference images, which are obtained by subtracting the results at 0 T from those at 11 T. 
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FIG. S6: (A) Difference between the conductance maps at B = 0 T and 11 T; Sg(r,E,B = 11 T) = g(r,E,B = 11 T) — 
g(r : E,B = 0 T). Black dashed circles mark the observed vortices. (B) Low-pass filtered image of (A); LPF-(5g. A cut-off 
wavevector is chosen to be 0.05 x 27r/ao in the filtering process. Corresponding spatial resolution is indicated as a black 
circle and an arrow. (C) Azimuthally averaged cross-correlations for pairs of [LPF-(5g, Ai] (Green curve), and [LPF-5g, 
Z{ r, e = 1, B = 0T)] (Red curve), and [LPF-£g, As(r)] (Blue curve). (D)-(F) Ai(r), Z( r, E, B = 0 T), and As(r) at B = 0 T. 
Dashed circles in these figures represent the locations of the observed vortices. 
























